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Abstract 

We (re)cosider the propagation of small disturbances (sound waves) in the presence of a pinned irrotational 
vortex in a superfluid with the help of the formalism of acoustic spacetimes. We give closed formulas for 
I the scattering angle for sound rays, formulate the sound-propagation problem in the Hamiltonian form, and 

o 

C\| discuss the form of boundary conditions at the core of the vortex, where the Hamiltonian has a singular 



point. The wave equation is simplified to a single ordinary differential equation of Mathieu type. We give an 



> 
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' extensive discussion of perturbations localized close to the core, which are similar to what is known as the 

. Kelvin waves. The spectra of modes depend strongly on the type of boundary condition employed close to 

,_i the vortex core. The existence of the gapless Kelvin mode, which is one of the modes with angular number 

^ -1, is usually discussed in the context of unpinned vortices in superfluid helium or rotating Bose-Einstein 

Q I condensates. We prove that this particular mode is absent if the vortex is pinned, and consequently one must 

d • discuss the full family of modes in this case. The question of whether or not the acoustic spacetime admits 

1 an ergoregion turns out to have a decisive (qualitative) influence on many aspects of sound-propagation 
i phenomena. 
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I. INTRODUCTION 



A. Background 

Vortices provide quantum fluids with a way to carry angular momentum jl, 2|. In the simplest case 
of a vortex line the stationary velocity field, v(x), is purely azimuthal and depends as |v| = const /r 
on the distance r from the line. This flow is locally irrotational but nonetheless possesses a finite 
global circulation F = J v ■ (ix. For fluids with phase coherence F must be an integer multiple of 
h/rriB, where itlb is the mass of the bosonic constituent of the fluid, e.g. the mass of a single helium-4 
atom for the superfluid helium-4. A way to probe rotating superfluids is provided by sound waves, 
jsj. In investigations of sound choices need to be made with regard to the type of description of the 
fluid and the structure of the vortex core (the innermost part of the vortex). Initial research on this 
iopic was mainly concerned with the superfluid helium-4 described by equations of fluid-dynamics 
3ue to the growing experimental accessibility of quantum vapors (BECs) the recent works 
tend to deal with this type of fluid, employing various approaches stemming from the Bogoliubov 
theory, js-^. The correspondence to earlier research is provided by (but also limited to the range 
of applicability of) the relation between these detailed descriptions and the fluid dynamics. 

In the case of experiments with vortices in helium the typical conflgurations consist of cylindrical 
cryostats fllled with the superfluid, with the possibility of placing a wire (or rod) on the axis of sym- 
metry of the cylinder. For not too large angular momenta of the fluid the velocity fleld corresponding 
to the stationary flow is exactly the one of the vortex line, while for higher angular momenta a vortex 
lattice is formed [3;] . The interesting problem of flnding the most stable stationary state of the fluid 
for given radii of the cylinders at a given teraperature and angular momentum of the fluid has been 
intensively investigated by many authors js-lOj. 

On the other hand, in the case of BECs, the typical configurations consist of trapped rotating 
clouds of supersaturated vapors. For repulsively interacting constituents the vapor settles down to a 
metastable state, which is well described as a weakly-interacting Bose gas [lllll2|. The type of the 
most stable configuration of a trapped gas with a given angular momentum depends strongly on the 
form of the trapping potential (the potential seen by the constituents in the plane perpendicular to 
the axis of rotation): for harmonic trapping an arrangement of single vortices (each carrying a single 
quantum of circulation) is the stable configurations for not too large angular momenta. Harmonic 
trapping has the special feature, that above certain angular momentum the BEC cloud begins to 
shed atoms radially. However, for trapping potentials "steeper" at large distances from the axis 
larger angular momenta can be accommodated by the cloud. When this is done, a configuration 
with a giant vortex in the middle of the trap is the expected final configuration at large angular 
momenta [9|]- A single vortex (simple or giant) provides the velocity field of a line vortex. 
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The axis of rotation of an idealized vortex line is a singular place because |v| = const /r cannot 
continue to r = 0. Dynamical models of vortex cores regularize this singularity. On the level of 
fluid dynamics one may distinguish two typical types of models: for superfiuids, if an equation of 
state is given (and the symmetry of the flow is assumed) then we are lead to a vortex model with 
an interior surface there, where the pressure reaches the lowest allowed value (the saturated vapor 
pressure). In this way an irrotational everywhere regular velocity field is obtained. Alternatively, 
for normal fluids, one assumes (as in the e.g. Rankine- 13|], or Lamb-Oseen 1^ models) that 
there is a distribution of vorticity mostly localized in the innermost part of the vortex. In this 
way the irrotational const /r profile is not continued towards the core, but rather the fluid velocity 
remains finite everywhere. Such vortex-core models are apparently adequate for the description 
of real phenomena in aerodynamics, jl^. On the other hand, in the case of rotating quantum 
vapors (BECs) there exist natural descriptions of the innermost region, alternative to fluid dynamics, 
ll2[, perhaps the simplest of which being the Gross-Pitaevskii description. While in regions 
where the condensate density is large and slowly varying in space the BEC can be described as 
a (super)fluid, this is not anymore the case in the core. Even though the equations of the fluid 
dynamics fail there, the mean-field Gross-Pitaevski (or other) description is still valid, and essential 
characteristics of excitations of the vortex can be computed reliably (numerically), see e.g. 
While these characteristics are very important for experimental verifications, they are nonetheless 
specific to the case of weakly interacting quantum vapors. In this paper we chose to work within 
the universahty class provided by equations of fluid dynamics. We show, that qualitative difference 
is expected for sound propagation phenomena once the region to the interior of Tc = T/cs (with Cg 
denoting the local speed of sound) is accessible to the fluid. 

In the variety of configurations outlined above the problem of sound propagation has been consid- 
ered by many authors. The simplest approach, also used in this work, is to employ the macroscopic, 
fluid-dynamic description. The problem is further simplified if the perturbation of density and ve- 
locity can be assumed small (linearization). Care must be exercised with this assumption if the 
position of the vortex is allowed to move as a result of the perturbation. If the velocity the flow 
to be perturbed is large close to the core (i.e. if the modification of const jr profile is localized at 
small r) then even a small shift of the axis of rotation may mean a large perturbation of the velocity 
field. In this work we consider only pinned vortices, for which the above obstruction does not occur. 
The other way to address the above problem is to consider regularized core models 14l-ll6| with 
sufficiently large cores, so that the background velocity is not too large anywhere. The study of 
perturbations naturally splits into the consideration of sound scattering and of perturbations bound 
to the vortex. 

For the scattering the natural problem is to derive the phase shifts in the partial-wave expansion 
of perturbations. Approximate approaches to the problem start with the work of Fetter j4|, and were 
continued with growing sophistication. Some renewed activity in the field appeared in the last decade 
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due to the connection to the lordanskii force problem, [is], 0, 13], and the controversy associated 
with it. We remark, that the general consensus on independence of scattering characteristics (for 
small frequencies) from the details of the vortex core can easily be shown unfounded in the case of 
a pinned vortex admitting an ergoregion. We do not focus on this problem, however, in the present 
work. 

As to the study of disturbances localized close to the vortex core it is important to note that 
generally - for each model of the vortex core - there exists a large family of such. The characteristics 
of the allowed modes, such as the dispersion relations for the (many) branches existin g fo r each 



angular number m, depend strongly on the type of regularization of the vortex core 14l4l6|. This 
family is usually referred to as the family of Kelvin waves. In many cases, there exists an important 
mode of vibrations with the angular number m = — 1, positive frequency u, and the gap less dispersion 
relation u{p) ~ p^ln(l/p) for small wavenumbers p (wavenumber along the vortex). This mode is 



usually referrec 
for a given p, 



to as the Kelvin mode. 
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19| . Modes with |m| = 1 (there are many such modes 
16|) generally put the cores of unpinned vortices into motion. This has lead 
many authors to seek for simple explanations of these modes on the phenomenological level (e.g. by 
considering forces acting on the core, as if it were a string with a specific tension, 0, Q). It should 
)e remarked, that a version of the Kelvin mode has also been apparently found in BEC systems 
sl. [lol. although it is puzzling why none of the other modes of the family (which is usually very rich) 
has ever been reported to be found in this context. 

Because for vortices pinned on sound-reflecting wires a qualitatively different dynamics in the 
innermost part is at work, and because such vortices do not allow for deformations of the core 
(vortex line cannot be considered as a deformable string), we shall refer to the modes localized close 
io the core/wire discussed in this paper as the quasi Kelvin waves (or whispering-gallery modes. 



20|). This is further justified by the fact that, as we shall show, the familiar Kelvin mode, with its 
characteristic gapless dispersion relation, is absent in the system for two types of one-parameter 
boundary conditions we consider. 

The main purpose of this paper is to report on new results for sound propagation in the presence of 
a pinned vortex line. We use the fluid-mechanical description and employ the formalism of acoustic 



spacetimes, which has been under intensive development recently, 2l| . In this way the mathematical 
structure of the problem receives additional attention. The sound fields are parameterized by a 
single scalar potential, and a simple Hamiltonian structure is available Q]. Quite remarkably, many 
technical aspects of the problem turn out to be related to other important areas of theoretical physics, 
including: stability of (classical and quantum) fields in the presence of ergoregions, behavior of fields 
in the presence of strongly binding potentials, and the Aharonov-Bohm scattering. While the last 
relation is well known among the authors Q]' point its limits once the acoustic spacetime 
contains an ergoregion. To the best of our knowledge the following results are new: exact formula 
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for the scattering angle in ray-acoustics approximation (section [IT]), discussion of the general form 
of boundary conditions at the core ("core conditions", section |V]), (semi-numerical) determination 
of spectra of bound states (quasi-Kelvin/whispering-gallery waves) localized close to the vortex for 
Neumann and "core" boundary conditions (section IVI A 2[) . and the discussion of various sections 
of the parameter space displaying the characteristics of these states (section IVI Bl) . The recurrent 
theme, throughout the paper, is the the qualitative change of the character of the problem once the 
acoustic spacetime contains ergoregions. Some of the problems related to stability of the flow are 
formulated, but left to a future investigation (see section IVI C|) . 

The paper should be regarded as an attempt to strengthen the understanding of the mathemati- 
cal structure of the sound-propagation problem for perturbations of the vortex flow. The additional 
features reported here in the fluid-dynamical approach (e.g. the specific characteristics of the full 
family of quasi-Kelvin waves) may or may not find their counterparts in the more sophisticated de- 
scriptions of concrete superfluids, such as the Bogoliubov descriptions of weakly-interacting rotating 

gJ . In section IVIII we present a short discussion of the possible exper- 



Bose-Einstein condensates 
imental regimes, where these features mig. 
is considered in a greater detail elsewhere 



it be detectible. The concrete case of superfluid helium-4 



20|. In our opinion, problems of superfluid rotation are of 



such importance that progress achieved even in very simplified models may turn out very valuable 
for future research. 



II. ACOUSTIC RAY DESCRIPTION 

In this paper we investigate propagation of sound in the presence of an irrotational flow of a 
vortex line. The vortex is assumed to be pinned (no deformations of the line allowed), and the 
density perturbation is assumed to be "small" (linear regime). The propagation of sound for such a 
flow can be derived from a single partial differential equation for a single scalar potential, ip. The 



velocity- and density-perturbations are expressed by, [22|, |23[ , 



6vi = dii), (1) 
6p = -^{dtiJ + v'd,ij), (2) 

where Cg denotes the local speed of sound, v the local velocity field, and p the local mass-density of 
the background fluid configuration. The potential ip fulfills the d'Alembert equation associated with 
the so-called acoustic metric, the form of which is 



ds^ = ^ [-c^dt^ + 6ij {dx' - v' dt) {dx^ - dt)] , (3) 
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In the stationary, cylindrically-symmetric case of an irrotational vortex line it reduces to 

ds^ = — [-fdt^ - 2k dtd^ + r^d^^ + dr^ + dz^] , (4) 

where f = cj. — K^/r^ and f^ = v,p = ^ = ^f vrfx is a constant of dimension [/t] = [cm^/s]. Positive 
K corresponds to a fluid rotating in the positive direction. This constant, together with a reference 
value of Cs, e.g. at oo, allows us to construct a length scale, Vc = n/cg, and a time scale tc = n/cj.. 
In what follows we will use dimensionless quantities t/tc and r/r^, denoted by the same letters {t,r) 
for brevity. From now on we will make the simplifying assumption, that Cg and p do not depend on 
the distance r from the axis. Whether or not this is a restrictive assumption depends on the type 
of superfluid considered. We discuss the deviations from this assumption for superfluid helium-4 
in chapter IVIII The case with variable p and Cs is considered in a separate work, [20]. With the 
assumption, and the rescaling, we are lead to considering physics of scalar fields in the spacetime 
with the metric 

ds^ = -{1 - ^)dt^ - 2dtdip + r'^dip^ + dr"^ + dz"^. (5) 

The ranges of coordinates are standard, save for r, which in some cases will be restricted from below 
by ro. 

The first task is to determine null geodesies, corresponding to the acoustic rays (ray acoustics 
approximation) . The equations of motion can be derived from the Lagrangian 

L = \ {-fP - 2t0 + rV' + r^ + z^) = l XaX^ = 0, (6) 

where / = 1 — ^2, and where differentiations are performed with respect the parameter r of the 
geodesies, as x"" = x"(r). The last equity is due to the null character of the rays. The spacetime 
contains an ergoregion, if points with r < 1 are accessible. This may, or may not, be the case, 
depending on the experimental configuration. If the vortex is pinned on a central wire, then the 
(rescaled) radius of the wire provides the lowest value of r, the tq. The problem of integrating the 
geodesic equations reduces to quadratures, because there are four first integrals: the momentum 
along the vortex, p = z, the "angular momentum" J, the "energy" E, 

J = r^^- i, (7) 
E = fi + if, (8) 

together with the Lagrangian L. Note that objects moving along the the zero-angular-momentum 
curves, J = 0, possess the (coordinate) angular velocity ^ = and thus remain at rest with 
respect to the fluid particles (i.e. co-rotate with them). 

With the help of E, J, and the Lagrangian we obtain the identity ^{—Ei + Jip + P + z^) = 0, 
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which is fulfilled along the geodesic. The reparametrization invariance of null geodesies allows us to 
set J or E to any value; we chose E = 1, leading to t = 1 at r = oo (provided, geodesies reach that 
point). The velocities (t, 0) are eliminated by 



t = r — J, 
. 1 + fJ 



t 



(9) 

(10) 



(The value of t could, as it stands, decrease [43[ along a null geodesic with J > if this geodesic, 
coming from r = oo, crossed the radius r^aus = This does not happen, however, as we shall see 
below.) 

The only equation to solve is the radial equation, which now assumes the form 



J 



{l + fJ)-J-p\ 



'111 



A. General properties of rays 



Let us consider rays which reach r = oo and call them scattering rays. The limit r — )■ oo of the 
radial equation and of Eq. ( !T0|) yields 



1 + J 



(12) 



Consequently, it is possible and advantageous to introduce the standard scattering parameters. All 
scattering rays are specified by: the projection of the velocity on the radial direction. 



def . / ^ 

f = Too = cos(aj 



and by the "impact parameter" , 6, defined in the way familiar from scattering theory 

hv 



(13) 



hv = l + J, 

(at r = oo). We obtain the following form of the classical scattering problem: 

{J+l)+v\ 



(14) 



(15) 



Let us mention, that the second type of rays appearing is obtained by taking \p\ > E = 1, that 
is with v"^ < 0. These rays do not reach r = oo and therefore are rays hound to the vortex. Setting 
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£ = v"^ (and allowing for negative £) we investigate the polynomial in ^ on the right hand side of 

r^ = J^(ir-J(J + 2)(i)2 + ^, (16) 

which is just a rearrangement of f[T^ . Zeros of this polynomial correspond to peri-/apocenters of 
the rays, and motion can only take place in regions where the polynomial is positive. For £^ > 0, 
motion will always be possible at the two limits: r — > oo and r — )■ 0. Setting ^ = w we note, that 
the product of the zeros is equal to £/J'^. For the bound rays, with £ < 0, exactly one zero of 
the polynomial is always positive, say (ui+), and it corresponds to exactly one positive root for K 
Therefore, the bound rays are coming from r = 0, reach the apocenter, and fall towards r = 0. The 
case of scattering rays, with £^ > 0, requires a closer look. 



B. Determination of orbits for scattering rays 

In this section we consider the scattering rays, and determine the form of their orbits, i.e. the 
functional dependence of 1/r = u on the angle ip. Quite generally we shall see, that geodesies 
starting from r = oo either reach a pericenter (and scatter back to oo) or propagate towards and end 
on the vortex core. The problem is integrable in terms of elliptic functions. The range of angular 
momentum J leading to absorption is, as expected, asymmetric due to the rotation of the vortex. 

We start with 

du dr dip /r T 1 T 2n /i 

r = — — ^ = -u'[J+l- Ju^] (17) 
d(p du dr 

and find from ( IT6l) 

{u'y [J + 1 - Ju']2 = [J + 1 - Ju^] [1 - Ju^] - ( J + 1) + v\ (18) 

Thus our problem is the classical Jacobi's inversion problem for u{ip) on a Riemann surface of genus 

du[{J+l) - Ju^ 



one, 

'•"0 \( T J_ _ T,,2 



with 

Pi{u) = - J{J + 2)u^ + v^. (20) 

In order to compute the scattering angle A we set Uq to the lowest positive root of P4 (pericenter), 
where we can assume = 0, and compute the integral up to m(A) = (corresponding to r — )• 00). 
The total scattering angle, 9, will be 6' = vr — 2A. Before establishing the dependence of A(f , h) let 
us first consider the case when no scattering takes place, because there is no pericenter, an the ray 
propagates to the vortex core. 
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Fall onto the core 



The polynomial P/^{u) is bi-quadratic. Let w = v?. Both roots of P2{w) = Pa{\^) have the same 
sign due to w+W- = v^/J"^. In two ranges of parameters (J, f) none of the roots, w±, is real and 
positive (i.e. they do not lead to real values of u, and therefore the ray continues r = 0): they are 
not real unless (J + 2)^ — 4f ^ > 0, and they are real but negative for J G [2f — 2, 0]. These ranges of 
parameters have been depicted in the Fig. [H For other values of (J, f) there are two real, positive 
roots, and thus four real roots of Pi{u), two of which are positive, say u±. The rays starting at 
r = oo get reflected at m_, while the ones starting at r = get reflected at m+. 




-5 -4 -3 -2 -1 1 -15 -10 -5 5 10 15 

J b 



FIG. 1: The range of parameters for capture of rays coming from oo: [Left:] {J,v) plane [Right:] (6, 
plane, with hv = J + 1. For small v, that is for rays with small velocities Too = v (propagating along the 
vortex line) the capture occurs for almost all values of the impact parameter h. 



By considering the limiting form of the radial and axial equations we find the following asymptotic 
behavior of rays falling onto r = 0: 

r'-^^. (21) 

from which it turns out, that the asymptotic trajectory is an (independent of J) hyperbohc spirah 

d/T 1 

- — =d!f ^!p = ifQ+-. (22) 
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FIG. 2: Exact (numerical) solutions for sound rays. [Left:] Acoustic rays coming from r = oo, which (1) 
get scattered by the vortex, (2) get absorbed by the vortex. [Right:] Rays emitted from the core which get 
reabsorbed. Paths have been plotted only up to the peri-/apocenter, where appropriate. Also visible is the 
boundary of the ergoregion, located at r = 1. The flow in the vortex is positive-oriented. 



Exact results for the scattering angle A 



Appearance of a quartic in the radical indicates a solution in terms of complete elliptic integrals. 



We recall the definitions, 



K{m) 



2i- 



dt 



lo V(l-t2)(l-mt2)' 
and note the following identities 



E(m) 



1 — mt^ 

l-t2 



dt. 



(23) 



dx 



y/x'^ - Ax^ + 1 Xh 



K{xt), 



Vx^ -Ax'^ + 1 



x+[K{xt)~E{xt)], (24) 



where A > 0, X- is the smallest positive root of the radicand (which has, per assumption, four real 
roots), and = 1. Taking v out of the radicand and introducing a new variable z = u^/\J\Jv 

we obtain expressions of the above form with the result: 



A(J,.) = ^^ir(/)-^^"('^)^~ r....4^ ...4 



^v\J\ z+ 



[K{zt)-E{ 



z_ 



(25) 



where z± are (positive) roots of 



A ± - 4 



with A 



J + 2 



sgn (J) > 0. 



(26) 
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The total scattering angle is just 6 = n — 2A. The scattering angles have an interesting expansion 
for large values of the impact parameter b = {J + l)/v: 



where the negative values of b correspond to rays coming from below the x axis in the left part of 
Fig. [21 The symmetry of the first two terms, as well as an overall focusing character of the vortex 
scattering [43], drew attention of authors before 25|. On the "positive absorption edge" (J = 0) we 
have A — ?■ +oo, while A — )■ — oo on the "negative absorption edge" (J = — 2f — 2). Note also, that 
due to f ^ < 1 there holds 

W- = <-j = W,aus, (27) 

and therefore geodesies with J > coming from oo, which potentially could enter the region with 
if: < 0, get refiected before reaching that region. 



III. MATHEMATICAL STRUCTURE OF WAVE-PROPAGATION PROBLEM 

In this chapter we consider the propagation of (density) waves in a fiuid exhibiting a line-vortex 



fiow following the acoustic spacetime approach 22l, |23[. In this approach one notices that, as long as 



the background fiow is irrotational, the perturbation {6p, 5v) is parameterized by a single potential 
potential tp. This potential fulfills a single partial differential equation having the form of the 
d'Alembert equation fl53p in the acoustic spacetime, the metric of which in our case is given by 
While the most straightforward approach would consist of writing and separating this equation in 
order to determine the available "sound modes" , we first present the mathematical structure behind 
the problem. On the one hand this provides us with connections to other problems of theoretical 
physics, and on the other helps to avoid missing important caveats which, in our opinion, have not 
received due attention in the literature. 

A. Geometrical preliminaries: constant t foliation 

The acoustic spacetime has the line element 

ds'^ = -f dt^ -2dtdip + r^ d^p^ + dr'^ + dz^ (28) 



11 



with f — 1 — r ^. We have: g =^det (gab) — —r^, and (displaying (t, ip) components only) 

We may distinguish t = const surfaces Ct, which foliate the acoustic spacetime. The unit normal 
vector N"- orthogonal to these surfaces, Na = —da{t) is always timelike. The Killing vector of 
^-translations — {dtY may be orthogonally decomposed as 

r« = aAT" + (30) 



with 



finally 



1 

a 



-^ = 1, = (0,-^,0,0), (^,-1,0,0), (31) 



iV" = (1,-^,0,0). (32) 



The vector fields: ef. = N"- and e| = — r^a = (0, ^, 0, 0) form an orthonormal tetrad. 
For every covariantly conserved current J" by Gauss theorem we have that 

g = / dS^Ja (33) 
Jct 

is independent of the foliation-parameter t, provided the fluxes through boundaries of Ct are either 
absent or integrate to zero. The volume element of Ct is dS"- = \/hd^y iV", which for the special case 
of the vortex spacetime is just d^X(5^, because the internal geometry of Ct is fiat: 

ds'^\c^ = habdx°'dx^ = r'^dip'^ + dr"^ + dz^. 

Note, that in our case the determinants fulfill —g = h = r"^. 



B. Sound propagation as a linear dynamical system 



1. Function spaces 

In this section we provide some general arguments intended to fix the mathematical context of the 
subsequent more concrete investigations. Let us consider generally the space of "initial data" of the 
evolution of the scalar potential '0(i, x) in time t. Because the equation governing the evolution (the 
d'Alembert equation) is of second order in 9t, the space of initial data on Ct will necessarily contain 
the information about the potential and about its "time-derivative" p = N"-da'4>. The standard 
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structure of a linear dynamical system associated with this problem will be sketched below. Firstly 
we introduce the space of real solutions of the wave-equation, S^. Because of the uniqueness of the 
time evolution, this space is equivalent to the space of initial data, i.e. to the space of elements of 
the form $ = {ipjp) on any of the Cauchy surfaces Ct. In what follows we will make no distinction 
in notation between the space of solutions and the space of initial data (phase space). We equip the 
structure with the Poisson bracket 

$2] = dS''U<^>,, $2] = j d'x (7A1P2 - Pi^2), (34) 

where the rightmost form is specific to the vortex acoustic spacetime, and ja stands for the covariantly 
conserved current 

Ja[$l, $2] = $l«9a$2 - 5a$l$2- (35) 

We also consider the space of complex solutions, or complex initial data, Sc, equipped with the same 
Poisson bracket (but with the first argument complex-conjugated). We mention, that the dynamics 
of the system will be governed by a Hamilotnian of equation ( l38i) . 



2. Relation to quantum fields in curved spacetimes 

The above structures are sufficient from the point of view of classical field theory. At this point it is 
worthwhile to consider also the deeper problem of quantum sound fields propagating on a given fiow 
of fiuid, such as the field of phonons in a BEC. There are reasons to believe that quantized density 
luctuations in fiuids with phase coherence do exist and behave as a quantum field (see, however, 
2!). Indeed, this type of reasoning provides the basis for sonic analogs of Hawking-, Unruh- and 



other phenomena traditionally restricted for quantum fields in fiat and curved spacetimes 21 1 



Quantization of fields usually consists of two steps: the relatively easy one 45|] concerned with 
the formulation of the algebra of field observables, and the more difficult one concerned with a 
construction of a well-behaved representation of this algebra on a Hilbert space (such as the Fock 
space). The construction of a ground-state representation for fields in stationary axisymmetric 



spacetimes without ergoregions was given by Kay 26|]. Technically what is needed, is a Schrodinger- 



type description of the evolution of the "positive-frequency" solutions of the field equation. In 
order to achieve this one forms a complex Hilbert space of (one-component) functions H, usually 
known under the name of one-particle Hilbert space. The elements of H are meant to correspond to 
the "positive- frequency" solutions in Sc, with different choices of "positive frequencies" leading to 
quantized-field representations related by a Bogoliubov transformation. The time-evolution should 
be represented in H as a unitary operator U 

In simple cases, such as the propagation of fields in Minkowski spacetime, the way to project 
elements of Sc into "H is provided by an operator j : Sc 1-L. For a solution "of frequency w" the 
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action of j would have the form 



(36) 



where p = dtip in that case. In our case, however, the frequency corresponding to an element of 
the initial data space iSr is not yet defined. The standard definition of the frequency is that it is 
the eigenvalue of the Hamilton-operator, H, which is t-independent for the stationary spacetimes. 
What is needed for the definition of the projection operator j, and thus also for the construction of 
the ground-state representation, is an operator h such that hip = ojip, hp = up for an eigenvector of 
the Hamiltonian H of the form {ip,p) = Then, the operator j will fulfill 



j[$(t)] =e^^*j[$(0)], i.e. U{t) = e 



iht 



(37) 



In the simplest case of a scalar field on the Minkowski space, the choice is h = V^, but we do not 
know if an appropriate choice can be made in the case of stationary spacetimes with ergoregions. We 
note, that the investigation of quantum sound fields can be carried further using methods developed 
by Kay 26|, although we stress, that these methods will not be available if an ergoregion is present. 
In this latter case no ground-state representation of the above type is known, and only more abstract 
solutions of the problem are available 
explored. 
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28|). the physical interpretation of which has yet to be 



3. Time evolution 

Returning to the classical fields, and denoting by T(t)$o the time-evolution of a phase-space 
vector $0 ^ Sc we derive the form of the Hamiltonian 



d_ 
di 



H(^Q = -a2A(^o, with if 



i=0 



i(5d i 
^iV^ i/3d. 



(38) 



where for technical reason we have distinguished the operator 



A 



(3% 1 



(39) 



while (72 stands for the Pauli matrix. Anticipating the separation of the angular dependence via e'"*'^ 
(with m G Z), leading to the partial-wave- expansion, we give the form of A for partial waves: 



72 im 



A 



(40) 
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where is the Laplace operator (of a flat 3-dimensional space). 

We note the connection to a technically very similar problem of complex scalar fields in electro- 
static potentials, as investigated in jg^l- Setting the (external) vector potential to zero (^ = 0) 
there, and denoting the mass of the field by /i, we obtain 

_J ^ I. (41) 



where V is the scalar potential. The correspondence between 27|] and our situation is therefore: 
V = ipd, which for partial waves reduces to = m/r"^. Thus, from the hamiltonian point of view, 
the behavior of (complex) sound potential in a vortex-background is equivalent to the behavior of 
a charged scalar field in strongly attractive/repulsive (depending on m) potentials. Such problems 



have been investigated in the literature before, see e.g. [29| and the appendix B of |28l|). 



For the linear space of Cauchy data one may introduce a natural, positive-definite, product 

($,$) = j d^y^/h^^^ = Jd^x^^^, (42) 

and also the product resulting from the Poisson bracket between and $, which (conforming to 
the literature) we shall call "the norm of , 

= ($, <^)jc = -($ , a2$) =tJ^d^yVh (# - p^) = ^] = J d'^ - ^) 1^-1', (43) 

(the last equity for ip of the form e~^'^'^~^^"^'^4'^{r)). The space Sq equipped with the product (., .) is 
a Hilbert space, while the space JC equipped with .)x: is a Krein (indefinite-product) space. This 



product corresponds exactly to equation (34) of |27| (with p = {dt — iAo)ip there), and to (2.15) of 



281]. The merit of introducing these additional function spaces lies in the special properties of the 



operators H and A. 



4- Properties of the Hamiltonian 

The Hamilton operator, H = — cr2^, which is explicitly given by 

(i(5d i\ f ^ i\ 

H = \ , for partial waves: H = \ ' (44) 

is a symmetric operator on the space /C, i.e. ($,i?\&)/(; = (if$, ^)x:. Additionally, in some cases, 
the operator A is a positive operator w.r.t. the product (., .) (that is ^ for all $ G Sc). 

In these cases one may proceed as Kay, and define the Hilbert (positive-product) space A with the 
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product {., .)a 

(vI/,<|.)^ = (M/,A<l>), (45) 

on which the Hamiltonian H = —a2A could be shown to be selfadjoint. In such cases, therefore, H 
has the standard properties of selfadjoint operators in Hilbert spaces, such as: the realty of eigen- 
values u, the orthogonality of the eigenvectors, and the completeness of the family of eigenvectors. 
All these properties, therefore, hold for sound modes propagating on a vortex background provided 
that the boundary conditions at the innermost part of the vortex lead to a self-adjoint H. 

In the cases where A is not positive definite, one is left with the Krein space /C, and with weaker 
results for H (see appendix |A]) , in particular: eigenvalues are not necessarily real (but do come in 
pairs of complex-conjugate numbers) and a broader class of vectors is necessary for completeness. 
As we shall see below, the distinction between cases - for stationary axisymmetric spacetimes - 
corresponds exactly to the distinction between spacetimes with/without ergoregions. 

C. Energy and angular momentum 

Let C,T = {dtY denote the Killing vector corresponding to time translation. The current associated 
with the notion of energy has the form 

= T\[ij]eT (46) 
where T"-^[ip\ stands for the energy-momentum tensor of the solution ?/'(t, x): 

Tab H'] = - \ 9ab d,^d,^, (47) 

and the index-brackets denote symmetrization, A(^ab) = |(^a6 + A^a)- The energy is obtained by 
integrating J!^ over any Cauchy surface Ct, 

E[ij] = [ dSaJ^ltP]. (48) 

JCt 

Note, that formally (if relevant boundary terms vanish) 

E[<t>^] = = -|(<l>^,a2iJ$^) = |(<l>^,iJ$Jyc = |c^||<^>.ir. (49) 

For scalar fields in Minkowski space the norm corresponding to positive frequency {u > 0) 

solutions is always positive; here this does not need to be the case (because positivity of -^[$0;] is 
not guaranteed). 
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We may also regard energy as a functional of the solutions of the wave equation ip{t, x), 



= \ y"(i3x{|V^|2 + |p|2 + ^[^a^p-pa^t^]} = i($,A$), (50) 

where the integration by parts, — J d^^'il){S/'^ip) = J d^x.\'Vip\'^, was performed. (All boundary 
conditions employed in what follows are such that the - a priori present - boundary terms do vanish.) 
The energy functional is evidently positive when r is bounded from below by a ro > 1, i.e. in the 
absence of ergoregions. 

The angular momentum associated with a field configuration ip, J[ip], will be found as the Cauchy- 
surface integral of the current associated with the Killing vector field corresponding to rotations, 

J[iP] = J d^xT^biV^e^ = J d^x Reipd^iP) =mJd'^{LJ- ^M', (51) 
with the last equity holding for solutions of the form e~^^~^^"^'^'ip{r), for which we also find 



J[ijj] = |m||*||l 



(52) 



IV. PARTIAL WAVE DECOMPOSITION 



The scalar d'Alembert equation reads 



-9 

1 1 

— - dr{r) dri> + g^^didjtjj — -dr {rdrip) + 



ijj^O (53) 



and is equivalent to a Schrodinger-type equation 



idt^ = H^, with = p) e K. 



(54) 



Thus, the problem of finding "modes" either simphfies solely to the eigenvalue problem for the 
Hamilton operator 

H^^ = 0;$,, (55) 
or (in case H is not a normal operator) leads also to equations for Jordan associated eigenvectors 



(56) 
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with a (potentially infinite) chain numbered by the index a G N and = The Jordan 

eigenvectors evolve in time via 

= e"^^* [$^^(0) - it^''-\0)] , (57) 

fulfill the equation fl5^ and therefore are associated with solutions of the homogeneous wave equation 
Hip = 0. If we had a finite-dimensional case at hand, the associated vectors would be necessary to 
obtain a complete set of vectors in the linear Hilbert space with the product (.,.). We presume this 
also to be the case in the present context, should Jordan vectors exist for the given H with chosen 
boundary conditions. 

Proceeding with our special problem, we separate the angular- and the 2;— dependence, as is 
usually done, by 

$(t, x) = e^(™'^+P")^(t, r), meZ, peR. (58) 

The domains of m and p are justified by single- valuedness of $ and self-adjointness of —idz, respec- 
tively. The problem is stationary, and for this reason we separate the time-dependence by 



where, however, w's are not a priori restricted to be real numbers. Let us consider the eigenvalue 
problem for H, eq. (!55i) . in a sector of fixed {m,p). The /C-product, eq. (H3l) . of two eigenvectors of 
H assumes the form 

($^1,<I'^J^= / rf^x(cJI + CJ2 - ^)^^a;2, (59) 



where <l>^ = {i^^, -i{uj - ^)V'a;)- 

The eigenvalue problem reduces to solving the radial equation for ip = ip^^j, which is the main 
equation of this paper: 

1 r , A/f2 rr)2 1 

(60) 







^ = 0, 




r 







with 



46| 



e = u''-p\ 



= m{m + 2u). 



Neither nor k"^ need to be positive; negative values of these parameters will correspond to 
physically distinguished situations. From now on we will use the notation 



Dt 



d? 1 d 
dr"^ r dr 



for the (often appearing) radial part of the two-dimensional Laplace operator. The frequency, u, 
as an eigenvalue of symmetric operator, H, on a Krein space does not need to be real. The radial 
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equation has the symmetry {m,u) — )■ (— m, — w), so that finding all positive frequency solutions 
would be enough. We note the singularity at r = 0, which in the usual case (i.e. if we were dealing 
with the usual Hilbert space of square-integrable functions on R^) leads to well-known ambiguities 
associated with selfadjoint extensions of the Hamilton operator 30| . 

If not for the scalar product, the problem would be of the Schrodinger form with 



m 

^4 



(61) 



V{r) 



with the potential V{r) diverging to — oo as r — )■ 0. 

Momenta along the vortex, p, are not necessarily restricted by p"^ ^ u"^. Such a restriction in 
usual scattering problems comes from the required behavior the innermost parts of the "space" , for 
example from a boundary condition at a fixed r, or from a regularity requirement at r = 0. We will 
investigate this problem in section |Vl 

For comparison with ray acoustics, we note the relations 



i = E- 



J 



(62) 



and 



{E'-p' 



J{J + 2E) , J 
r 



+ -' 



so that for apt analogy with ( !6T1) we should make identifications m i-^ J, u ^ E, and k'^ = E'^ — p'^ 
We also note the correspondence between ipip = {u — ^)ip and t of the acoustic rays. 



V. BOUNDARY CONDITION AT/NEAR THE VORTEX CORE 

The Hamiltonian is defined by its formal expression, ( H4|) . together with appropriate boundary 
conditions. We recall, that in the case of operators in Hilbert spaces choosing a "wrong" boundary 
conditions can make a symmetric operator (fulfilling {ip,Hx) = {Hip,x)) non-selfadjoint, opening- 
up the possibility of complex frequencies u. For a one- dimensional differential operator posing 
a boundary condition at a regular point of the corresponding ODE is relatively unproblematic. 
Moreover, the case of boundary conditions supposed to be posed at singular points of the ODE 



there is also a well established theory due to Weyl 3l|- In the vortex case, assuming the whole space 



r G [0, oo) is available, we have the problem of posing the boundary condition at r = 0. This point 
is a singular point of the radial equation, and the Hamiltonian is not given as a symmetric operator 
on a Hilbert space. 

Let us recall an important example. For the Schrodinger problem (in Hilbert space) with the 
potential V = —'^ Case proved jsO], that for small r the eigenfunctions necessarily have the asymp- 
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totic form ip ^ rcos + with the same value of the phase B for all to. Posing a boundary 
condition at r = is equivalent to a choice of B, which then becomes the essential property of 
the operator. This inhomogeneous condition replaces the more commonly employed homogeneous 
boundary conditions (such as the Dirichlet/Neumann conditions) at r = 0, and it is necessary for the 
mutual orthogonality w.r.t. (., .) of eigenf unctions to different eigenvalues (a "must" for symmetric 



operators on a Hilbert and Krein spaces) 47|. A reasoning, similar to what has been done in 30|], 
can be repeated in our case. 

Close to r = the radial equation assumes the form, which it always has for k = (i.e. in the 
case of waves the z momentum equal to the frequency, p = G M). Introducing m = 1/r we find, 

u^i;" + uip' + {m\^ - M^)ip = 0, (63) 



and thus 

ip = BesselA/(|^|/T), for r ^ 1 or for k = 0, (64) 

where BesselAf(x) stands for an arbitrary regular Bessel function {J,Y or H^''^) of order M. We 
recall that even for w G M the parameter M can be real or imaginary (the latter case occurs for 
m G {—2uj, 0), which is the wave analog of ray-capturej^sl). 

The asymptotic expansions of the radial function at r <^ 1 for the distinguished Bessel functions 
read 

Jm(^)^ cos(f + M^-^).y|, (65) 
YmO-^) ^ -sin (f + - ^) ■ . (66) 

For any solution of the radial equation fl60|) . ip, the asymptotic expansion at small r can be charac- 
terized by a phase B: 

^|J^cosB. Jm(^) - sinB ■ Ym^-^) ^ cos (f + - ^ _ 5) . (67) 

or (shifting 5 by f + ^): 

^^cos(M + 5).y^. (68) 

We now repeat Case's reasoning in order to reveal the type of boundary condition at r = necessary 
in our case. Let $1 and $2 be two such solutions corresponding to the frequencies Ui, U2, and 
characterized by phases Bi and B2 close to the core. Assume, that both correspond to the same 
values of the separation constants {m^p). In order to simplify matters at r = 00 we assume both of 
the functions to fulfill a Dirichlet boundary condition at some large r = R, and investigate them on 
the interval r G [0,R]. Multiplying the radial equations ( !60!l fulfilled by each function by the other 
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function, subtracting them, and integrating the result from to i? we obtain 



R 



{U2 - Wi)($i,$2)x;- 



(69) 







(More precisely, there should stand only the radial part of the /C-product on the r.h.s..) Because 
eigenfunctions /C-selfadjoint operator to different eigenvalues (ci;i,ci;2) must be /C-orthogonal (see 
appendix [X]) , using the asymptotic expansion at r = 0, eq. (1^ . we arrive at 



where on the l.h.s. a (regular) limit r — )■ has been performed. Thus, for each partial-wave problem 
(fixed {m,p)), the phases B must be independent of the frequency u. This is so unless we are 
dealing with a pair corresponding to mutually conjugate, complex frequencies, which - however - 
might appear in the problem. The boundary condition, as far as it can be fixed by the orthogonality- 
argument above, amounts to setting the phase B{m,p) (a function of m and p) to be independent 
of u. In real physical situations the values of the phases B{m,p) will be related to the physics of 
the innermost parts of the vortex, and will depend on the model of the core used, e.g. on the type 
of the fluid in motion. In the subsequent chapters we will illustrate the fact, that the choice of these 
phases does influence essential physical characteristics of sound-propagation on the vortex. We will 
do so, by taking the B{m,p) to be independent of {m,p), essentially choosing either the Bessel- J or 
Y asymptotics for all m and p. We do not expect this type of choice to be related to any particularly 
distinguished physical situation. 

VI. SOLUTIONS TO THE RADIAL EQUATION 
A. Structural results 

1. General considerations 

A general-enough consideration of the radial equation, eq. fl6Up . assumes only m G Z and p G M. 
It is important to restrict the number of parameters to an essential set. Recalling, that a priori G C, 
and thus k"^ = uj"^ — p^ is in general a complex number, we introduce a new variable = l^l^"^, 
X G [0, 00) and transform the radial equation (160|) into the form 



where = m(m + 2u)), and \mk\ = K^. We note, that by a simple transformation of variables the 
equation can be shown to be equivalent to the Mathieu equation, taken along various contours in the 



m\ sin(5i - B2) = {co2 - Wi)($i, $2)^ = 0, 



(70) 



(71) 
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complex plane (see appendix [U]) . In the case of real frequencies the equation assumes the form(s 

1 



±1 + 



^ = 0, 



(72) 



where the + sign corresponds to solutions having an interpretation of scattering states, occurring 
for k"^ = u"^ — p"^ > 0, while the — sign to "bound states", possibly occurring for k"^ < 0. Generally, 
unless k"^ is real positive, exactly one of two linearly independent solutions will be exponentially 
damped at a; — t- oo. 

Let us consider the scattering, k^ > 0, case first. On the two extremes of the domain of x (0 
and oo) we drop the fastest decaying terms: or x"^ = u'"^, and obtain two Bessel equations with 
equal indices 



M2 



^ = 0, 
^ = 0, 



X > 1, 
u > 1, 



where u = x ^. 



2. Bound states in the (Er^,M^) plane 



Having established the general form of the radial equation we now study its distinguished solutions 



more closely. As remarked before, in the case of negative k"^ 

1 



-1 + 



x'' 



ij^ — p^, the equation has the form 



= 0, (73) 



with X G (0, oo). The solutions of this equation behave as modified Bessel functions {Km{Kx) and 
Im{Kx)) at X — 7- oo. An important question arises as to whether for a given boundary condition at a 
finite (or zero) x there exist values of M and for which the solution would be "bound-state-like" 
with the asymptotics containing only the exponentially damped Bessel function. One might expect 
to be able to answer this question with the help of the specific investigations of the Mathieu equation 



le.g. 
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33|). So far, however, we have not been able to locate the answer in the literature. 



For that reason we present here numerical results, from which important physical conclusions can 
be drawn. We shall use two types of boundary conditions: 

I) The family of "core" -boundary conditions provided by the orthogonality argument of section 
IVl The values of the constant B corresponding to asymptotic functions of the first (J) and 
second (Y) kind will be considered. We will only consider the (illustrative) cases, where the 
phase B is independent of m and p; otherwise the boundary condition would depend on the 
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position in the {K'^,M'^) plane. 

II) Recalling that in the acoustic spacetime formalism the perturbation of the velocity field is just 
the gradient of ip, we consider the family of Neumann boundary conditions posed at a finite 
r = vq (corresponding to a flow around, say, a central wire of radius tq). This allows us to 
consider the effect of the presence/absence of an ergoregion, which is included in the spacetime 
if points with r < 1 belong to it. 



3. Scaling transformations, and positions of bound states 

Both types of boundary conditions, the Neumann condition and the "core" condition f l68|) . can 
(for given values of regular parameters (m,p, u, tq or B)) be written purely in terms of the essential 
parameters/variables {x,M,K). For the Neumann condition we require 



dip 
dx 



(74) 



X=X(i 



(with the xq computable from ro, and from the other regular parameters), whereas for the "core" 
condition we require 

^ ^ cos(f + S)v^. (75) 

With this reformulation performed, we consider the radial (Mathieu) problem in the space of essential 
parameters. The results (such as positions of bound states in the (-ft'^, M^) plane) can always be 
recalculated to the space of regular parameters. An important ambiguity, however, occurs: because 
the number of essential parameters: (xq, M) for the Neumann condition, and (S, K, M) for the 
"core" condition, is smaller than the number of regular parameters, (rg, m,p, u) or (i?, m,p, w), there 
exists a one-parameter family of transformations for the regular parameters, leaving the essential 
parameters invariant. Recalling the relations 

a;2 = |A|r2^ K'^ = \mkl = m{m + 2u) (76) 

we see immediately that essential parameters remain invariant if (ro,mo, ^oji^o) is transformed into 

rx = Ari, mx = Ami, kx = jh, Qx = ^^i, (77) 

where 

n^ = l + ?!£x_ (78) 
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The only restriction on A is, that m\ should also be an integer. The frequency transforms according 
to 

= ^(i - A) + (79) 



4-. Positive /negative frequency hound states 

The A-transformations allow to flip the sign of the frequency of a bound-state solution. Intro- 
ducing an abbreviation /ii = ^ we consider 

^A = ^ = /^i(i-A) + i. (80) 

An inspection of the above function shows, that unless /ii G [—1, 0] (i.e. unless the (cui, mi) parame- 
ters fulfill mi G [— 2(X'i, 0], which is equivalent to ^ 0) the range of the function is the whole real 
set R, and therefore the sign of the frequency will be flipped there, where Ci;;,^ < 0. The change of 
the sign of oj^^ occurs for large A if yUi > and for small A if yUi < — 1. Thus a solution corresponding 
to cJi > 0,mi > can always be A-transformed to produce a solution with uj\ < Q,m\ > 0. This 
is true in the "existential sense", i.e. it allows us to state that such solutions would exist in certain 
boundary-value problems. The A-transformation for the Neumann boundary condition changes the 
radius tq (where the condition is posed), and therefore changes the problem considered. Solutions 
with the same sign of u and m rotate in the positive direction (co- rotate with the vortex). 

5. Negative norm states for positive frequencies 

For principal reasons in is important to investigate, whether positive frequency solution could in 
some cases correspond to negative "norms" : 

||$^||2 = 2 j d^^{u - (81) 

With this regard let us stress, that for positive the continuous nature of the positions of bound 
states in the (-ft'^, M^) plane allows us to investigate cases where = m^, with m G Z. Because of 
= m(m + 2uj) we may regard a; to be a linear function of indexed by the branch-numbering, 
discrete parameter m, 

uj^{M') = \ (Ml - rn) . (82) 

Then, a;m(M^) crosses zero at = m^, as we move along the bound-state-curve in (iT^, M^) plane. 
It is clear from fIST]) . that in regions of negligible u and positive m the "norm" of the corresponding 
solution will be negative. Therefore, states with slightly greater than the branch-defining m^ 
will have positive frequency Um^M"^) and negative norm | |$^^(a/2) | p. (Whether or not the norm of 
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such a solution will turn positive as we continue along the bound-state-line in the {K^, M^) is not 
clear.) 

On the other hand, the states with negative always correspond non-zero frequencies; the 
lowest possible frequency UmiM'^) for a given negative is u = \/—M^ (this happens only in 
case m = —u is an integer). The question whether or not negative norm solutions for positive 
frequencies exist for < does not appear to be answerable by general arguments of the above 
type. 

6. Bound states and scaling transformations 

The A-transformations allow us to tell, whether bound states are necessarily related to the presence 
of an ergoregion in the acoustic spacetime. We will consider the Neumann conditions, posed at tq, 
with ergoregions disappearing from the spacetime once the condition is posed at > 1. 

Let us take one solution related to a bound state, with the Neumann condition at ri, and consider 
the A-transformations. It is clear, that following the A-trajectory, the new position of the boundary, 
rx = Ari will - for sufficiently large A - be greater than 1. Focusing on a fixed "original problem" with 
the parameters (ro,mo, ko,uo) we notice, that for very large A's the following limits arise: mx ^ 1, 
A;^ ^ 1 that is p\ ~ ul, with the frequency ux = (and therefore |u;a| ^ 1). Such states, 

therefore, correspond to disturbances rotating in negative direction (against the vortex-rotation), 
with large angular quantum numbers, large frequencies and equally large momenta along the vortex. 




FIG. 3: Positions of the bound states in the {K'^,M'^) plane for core boundary conditions, together with 
distinguished lines M S N. [Left:] for the boundary condition ip ^ Jm{K/x), [Right:] for tp « Ym{K/x). 
The lines are zeros of the Wronskian of the exact (numerical) solution with the Bessel function Km{Kx), 
which is the asymptotic form of the exponentially decaying solutions at x = oo. We note, that none of the 
branches seems to enter/exist in the ^ quadrant, which suggests, that there are no bound states there 
(for the "core" boundary conditions). 

Let us remark on how figures such as FigJH are to be interpreted. Firstly, the value of Xq is the 



25 



0.05 0.10 0.15 0.20 0.25 



50 100 150 



50 100 150 

K'- 



FIG. 4: Positions of the bound states together with distinguished hnes M G N for the Neumann boundary 
condition at xq = {0.1, 0.5, 0.7}. Ah these plots are inequivalent in the sense, that they cannot be obtained 
from each other by a A— transformation. 

specific characteristic distinguishing between the spectra (lines in the (-ft'^, M^) plane corresponding 
to bound states). Secondly, the k = a/p^"^-u^ will be known, once Xq, tq and m are known, 

\k\ = \mr^r, (83) 

which allows us to write 

= m\^f. (84) 

This equation allows for multiple interpretations of a given point (fr^,M^) on the "bound-state- 
line", depending on the chosen value of m G Z. Thus, from given (i^^,M^,m) we compute the 
radius tq, and the frequency cOmiM"^): 

Furthermore, once Tq is known, \k\, and (together with the known Um) \p\ = ^/k"^ + can be 
computed. 

Points lying near = n"^ {n ^ N) are distinguished, because they can correspond to zero- 
frequency modes. For instance, from the Fig. HI with xq = 0.1 we find a bound state at = 1, 
^ 0.034, corresponding the zero-frequency, cJi(l) = 0, solution with |m| = 1 and 

ro ^ 0.54, p = ±{^f ^ ±0.034. (85) 

Thus, in this case, it is a problem with an ergoregion, and the (essentially static) perturbation has a 
form of a helical structure, with weak dependence, localized close to the vortex [ip ~ Ki{pr) for 



\m\xo 




— m] 
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large r). By considering states close to M^, with the same \m\, we obtain problems (characterized 
by '"o's) with similar solutions, but either co- (for > 1) or counter-rotating (for < 1) with 
the vortex. 

On the other hand we could have interpreted the above values of {K"^, M^) as coming from the 
second branch, \m\ = 2, for which we would obtain 



This solution, therefore, corresponds to a problem without an ergoregion, and the perturbation 
(localized very close to the vortex) is counter-rotating (sgn(a;m) < 0) and moving along the vortex. 

B. Location of bound states in the space of physical parameters; dispersion relations 

The scaling transformations in the previous section allowed us to formulate "existential" state- 
ments to the effect that the existence of one bound state for certain values of physical parameters 
To, u, m, p was equivalent to the existence of bound states for other values of these parameters. From 
the physical point of view it is, however, of interest to determine also the full "spectrum" of solutions 
for a given boundary problem. In other words: for a given type of boundary condition at r = 0, 
or at a finite r = tq, we shall compute the lines corresponding to bound states in cuts of the the 
{p,u,m) parameter space. This is a rather large space of parameters, and below we will attempt 
to extract the essential characteristics of the solution of the problem, that is essential features of 
the three-dimensional surfaces embedded in the four- dimensional space (ro,p, u, m). We restrict this 
numerical consideration to the boundary conditions of Neumann type posed at a finite r = tq. Only 
the case a; > will be presented, as the case a; < follows from the symmetry {u,m) o (— — m) 
of the problem. The form of the dispersion relations is qualitatively different for co- and counter- 
rotating solutions, as defined by the sign of m. (There are no bound states with m = 0, as the 
radial functions would need to be Ko{kr), k = a/p^ — w^, which however is a monotonic function 
and cannot fulfill the Neumann boundary condition at any tq.) 

For the counter-rotating solutions (m < 0) we obtain a number of branches concentrated 
just below the diagonal u = p (see Fig. E]). Solutions generally exist for all values of tq, but small 
frequencies are only possible in the presence of an ergoregion, i.e. if tq < 1. The group velocities 
approach (but are always less than) 1, meaning that packets constructed out of the wavefunctions 
corresponding to bound states would move along the vortex with almost the speed of sound. Given 
the observation, that each branch of solutions appears to begin at u; = p we may ask: what is the 
lowest frequency w of a solution(s) corresponding to various m for a given rg. In Fig. [6] we have 
plotted the positions of the branches of bound states in the {u, m) plane (with p = a; + 0.002). There 



ro ^ 1.09, 





(86) 
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FIG. 5: [Left:] Typical "band-structure" for bound states localized close to the core, propagating along the 
vortex with a momentum p, and counter-rotating with it for (ro,m) = (1.09, —2). [Right:] Rotated by 7r/2 
and extended view for the same values of (ro,?n.). 



is an apparent qualitative change between ro < 1 (ergoregion) and tq > 1 (no ergoregion): in the 
latter case no solutions with w ~ are possible. In the cases barely admitting ergoregions (rg ^ 1) 
the only branch approaching a; = is the leftmost one, and it does so for rather large |m| (see the 
central plot of Fig. |6]). For ro > 1 the solutions with small |m| are the ones leading to lowest w's. 



although not necessarily the m = — 1 mode. 
Kelvin mode in the unpinned situation. 



hese modes are all gapped, in contrast to the gapless 

3- 






FIG. 6: Positions of counter-rotating bound states withp = a;+0.002 for ro = {0.5, 0.9, 1.2}. Zero-frequency 
states exist only for ro < 1 and are related to leftmost branches in the plots. From ro > 1 up to ro ~ 8 the 
state with m = — 1 is the lowest-w state, while for larger ro this role is played by states with ever lower m. 



The co-rotating solutions (m > 0) exist only in the presence of an ergoregion. The branches 
are characterized by negative group velocities, all of which start at the diagonal, u = p, and end on 
the line a; = (see Fig. [7]). As the Neumann radius ro is increased (for a given m) the branches 
move towards small p, and some of them disappear. For each m there exist a critical radius ro, above 
which no co-rotating solutions exist. In the Fig. [9] we have depicted the positions of bound states in 
the (ro, m) plane, assuming fixed, small, values of frequencies and momenta, with u ^ p. Given the 



28 




5 10 15 20 25 30 5 10 15 20 25 30 5 10 15 20 25 30 

P P P 



FIG. 7: Typical "band-structure" for bound states localized close to the core, propagating along the vortex 
with a momentum p, and co-rotating with it for [Left:] (ro,m) = (0.1,5), [Central:] (ro,m) = (0.3,5), 
[Right:] (ro,m) = (0.1, 1). The plots correspond to a spacetime with an ergoregion; otherwise co-rotating 
bound states do not seem to exist. 

way position of bound states moves as ro is increased (Fig. [7]) an appearance of state at p ~ w ~ 
signifies the entrance of a new branch into the {p,u)) plot (dispersion relation). The rightmost line 
in Fig. [8] is the line of occurrence of the first co-rotating bound state, i.e. at a fixed m no bound 
states exist if tq is greater than the position of this line. For tq just smaller than 1, the co-rotating 
states exist only with very large angular numbers m. 
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FIG. 8: Position of co-rotating bound states with small uj and p [p = 0.11, u = 0.1) in the presence of 
an ergoregion as a function of the boundary radius ro and m. No bound states exist for tq > 1 (without 
ergoregion), and these that exist for large ro < 1 are necessarily associated with large m. 

C. Stability of bound-state modes (complex frequencies), relation to Friedman 
(ergoregion-) instability 

In the fluid-dynamical context, by considering small (test) perturbations of a given flow one can 
distinguish two types of instabilities. The first type develops, when the global energy of the fluid 
is lowered upon introduction of some perturbation. The second type comes about when modes of 
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perturbation with complex eigenvalues are allowed, that is fulfill the appropriate wave equation and 
the boundary conditions. An example of the the first, milder, type of instability is provided by 
the creation of excitations in superfiuid fiows exceeding the Landau velocity 3J]. It is natural to 
expect this type of instability in the presence of ergoregions because of the non-positivity of the 
energy density of perturbations there. However, it is difficult to estimate how rapid the dynamics 
associated with the instability will be, as the rate of creation of (negative-energy) excitations depends 
on external factors such as the surface roughness for capillary superfiuid fiows. In the vortex case, 
one can easily think of situations, where the fiuid will not encounter any factors efficiently generating 
negative-energy excitations. Invoking the spacetime analogy: it is difficult to estimate a priori the 
probabilities of Penrose- type processes, [35 1. 

Due_to the relation f H9|) negative energies are associated with positive-frequency, negative-norm 
statesj49|. Naturally, this cannot happen when m < (for counter-rotating modes), because the 
norm is explicitly positive definite in this case. For co-rotating (m > 0) modes we can generally assess 
that at least the energies of the modes close to a; = will be negative. This is always possible (see 
Fig. [7j) whenever such modes exist, but as we have argued it requires the presence of an ergoregion 
(see e.g. Fig. |8]). 

With regard to this first-type (negative-energy) of instability let us note the connection to the 
work of Friedman, where it is argued, that astrophysical spacetimes with ergoregions but without 



horizons are unstable [36|. The theorem proved there states, that unless a test, massless field can 
settle down to a non-radiative time-dependent state, the energy of the system in the bounded 
vicinity of the object will tend to — oo, radiating an infinite amount of energy to null infinity, 
provided the system starts from a negative energy state. The vortex spacetime provides a model of 
a spacetime where an analogous statement should also holdjs^. Note however, that our co-rotating 
bound states provide exactly the type of solutions the the "unless" clause would exclude: they are 
time-dependent, possibly negative-energy, non-radiative solutions of the d'Alembert equation in a 
stationary axisymmetric spacetime with an ergoregion. We do not see how, in the astrophysical 
context, the existence of the bound-state-like states can be excluded, but such states, decaying 
exponentially at large r (as in our vortex case) would violate the assumed asymptotic form of the 



solutions of the wave equation (Eq. (23)) in [2 



. Moreover, as the existence of bound states with 



negative energy in our work is seen to be related to the presence of the ergoregion, we doubt, that 



one can argue against the existence of such states 



in any context) by examining fall-off conditions 



of the field at null infinity alone (as is done in 36|]). If the completeness property for eigenmodes 



of the Hamiltonian turns out to hold in the ergoregion case, where it is a symmetric operator on a 
Krein space only, one will argue, that every wavepacket can be expanded into the eigenmodes, and 
the part associated with scattering states will be radiated away, carrying positive energy, while the 
part associated with bound-states will remain close to the innermost part of the system forever. 
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FIG. 9: Penrose diagram relevant to the theorem of Friedman, 36| . Part of the energy of a field configuration 
set up on the Cauchy surface Sq is radiated away through R, which belongs to the null infinity. Because 
the radiated energy is necessarily positive, the energy of the same field configuration on the (non-Cauchy) 
surface Si is lower. 

Let us finally comment on the second type of instability, associated with complex frequencies, 
which is typical of unstable configurations of fluids. By the known properties of symmetric operators 
in Krein spaces (see appendix \^ complex isolated eigenvalues of the Hamiltonian, (155]) . always 
appear in pairs of mutually conjugate numbers. The radial part of the Hamiltonian and the boundary 
condition associated with it depend, however, on three free, real parameters: p, vq or m. Let us 
call such an operator H^, where for the purpose of the argument we investigate the dependence of 
the spectrum on one of the parameters, say ^, with the other fixed |5l[|. Relying on the analogy 
with similar (Hamiltonian, not necessarily self- adjoint-) problems 27, 2^ S^], and on the experience 



with finite-dimensional operators, 38|], we expect that one signature of an appearance of complex 



frequencies would be if a pair of real frequencies (two lines in the (^, u) plane, say for ^ > ^*) were to 



"merge" at a point (^*,a;*), where the tangent line would be vertical, ^ = oo. Beyond this point 
(i.e. for ^ < ^*) only two mutually conjugate complex frequencies would remain. In the Fig. [TO] we 
present the dependence of u on various choices of ^ in what we observe to be the general behavior: 
no points of vertical tangent are to be found in the spectrum. We therefore conjecture that no 
complex frequency modes ever develop from a merger of real frequencies (in the vortex problem). 
However, the other way for complex frequencies to enter the spectrum is for them to enter from the 
continuous spectrum, i.e. from the region > \p\. Such an entry would leave no trace on plots of 
the positions of bound states in {C,,uj) planes (because on these plots a; G M per assumption), and 
requires further investigation. 

In order to give an example of how complex frequencies can manifest themselves in problems 
involving sound propagating on rotating background flows we recall the result of Sipp and Jacquin 
{igI - In this work an inviscid fluid forming a Lamb-Oseen vortex, T{r) = ^1 — e~^^^^o^, is 
investigated and it is found that many of the Kelvin waves are damped, i.e. possess frequencies 
with negative imaginary parts. Moreover, as one moves along a single branch (mode) in the {p, u) 
plane with an (initially) real frequency u, a non-zero damping rate suddenly emerges at some critical 
value of p. The relation of this result to our above argument is however not straightforward, because 
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FIG. 10: Typical positions of bound states: [Left:] in the (ro,w) plane (Neumann problem) for fixed 
{m,p) = (1,6) [Center:] in the (171,00) plane for fixed (ro,p) = (0.9,6) [Right:] in the {p,uj) plane for fixed 
(ro,m) = (0.2,-5). 



Lamb-Oseen vortex involves vorticity of the background flow and cannot be handled by the simple 
acoustic-spacetime description 52j. This description leads to the structure with H being a symmetric 
operator on a Krein space (see chapter [III Bl) . 



VII. ROTATING SUPERFLUID FLOWS ACCESSIBLE TO EXPERIMENTS 

Here we will recall the commonly available experimental setups, and judge their relation to the 
problem investigated in this paper. Attention will be paid to supersaturated quantum vapors (BECs) 
and to the superfluid helium-4. 



A. Supersaturated quantum vapors 

Scattering length, mass, interaction strength, vortex circulation for ^'^ Rh 

Common experiments with rotating BECs (see e.g. jsl, [lol) use doubly spin-polarized ®^Rb atoms. 
The relevant scattering length in this case is the the triplet one 53| at = 105 ± 4ao, that is at = 



5 • 10 ^cm. Using the atomic mass of 5^1 ^^Rb, tub = 87u = 1.4 ■ 10 ^^g we estimate the interaction 
strength 

, = ^ = 4.5 . 10-^« ^ (87) 

The quantum of circulation is Fq = / Vidx^ = 2TiVtp = ^ = 4.5 ■ 10~^ In the acoustic metric 
there appears 

a cm^ 

K = v^ = 7-10"^ . 
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(Some) Experimental values ofn, Cg, ^ and r^. 

Typical contemporary experiments involve condensate clouds with = 3 ■ 10^ atoms in a cigar- 
shape of characteristic dimensions 10^ x 100 /im^ = 10~*cm^ corresponding to a typical number 
density of n = 3 • lO^^cm"^. The typical speeds of sound are therefore of the order 



rriB 



For vortices we get the critical radius (radius of the "ergoregion" ~ 



K 



— = 7- lO^^cm = 0.7/im, (89) 



which is exactly the coherence length (below which the density and begin to vary with r) multiplied 
with 

^ = ^/V^2m^ = rc/y2 = 0.5/im. (90) 



The hydrodynamic description of BECs becomes inadequate on distances comparable with ^, and 
therefore even if the the core of the vortex could be approximately described by boundary condition 
of one of the types considered in this work, it would still correspond to an acoustic spacetimes 
without ergoregions. As a consequence of the relation = V2^, conceivable modifications of setups 
with a BEC- vortices (involving changes to number density n, interaction strength g, mass/type of 
the boson m^) would not change this state of affaris. On the other hand, the issue is open for 
pinned BEC-vortices or vortices in condensates trapped by anharmonic external fields due to the 
modification of vortex structure in such setups, [sl, 0]. Namely, e.g. for condensates trapped by 
potential growing quicker than (or for condensates trapped by cylinders with solid walls) a new 
picture arises The regular Tkachenko lattice of vortices begins to get replaced in the center by 
an empty core (giant vortex), where condensate density is zero. The velocity circulation around 
this core is finite, but the velocity of sound becomes negligible due to the vanishing of condensate 
density. Even though there are (usually) still vortices in the bulk, and the density is far from being 



constant everywhere 55|], the system exhibits characteristics of a spacetime with an ergoregion. 



B. Superfluid helium-4 

Single vortices 

The typical speed of (first) sound in the superfluid helium-4 is Cg = 2.4 x 10^^. Together with 
the mass of "^He atom, mne ~ 4m = 6.64 x lO^^'^g, this corresponds to k = = 1.588 x 10"'' — , 
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where Tq = is the quantum of circulation; with these values we obtain 

= 0.66 X 10~^cm (91) 

which is less than the van der Waals radius of the Helium atom, rvdw ~ 1-4 ■ 10~®cm. Consequently, 
single vortices in superfluid helium-4, when pinned on a wire to restrict their position, provide an 
arena relevant to this work, but with the Neumann (divided by Vc) greater than 1 (no ergoregion). 



Structure of a multiply charged vortex line; ratio v/c 

One would hope to create a model of ergoregion just by increasing the rotation of a vortex pinned 
on a wire. Let us contemplate this possibility for helium-4. General features of the structure of 
vortex lines can be inferred from the empirical, T = 0, equation of state (EOS) published long time 



ago by Brooks and Donnelly, 39|: 



3 

p = Y,Mp-por, (92) 

n=l 

with Aq = 560, Ai = 1.097 x 10^, A2 = 7.33 x 10"^, (with units giving pressure in atm), where po = 
0.145 g/cm'^ is the density at the saturated vapor pressure (SVP). The equation gives approximately 
correct compressibilities as well as the Griineisen constant Uc = ^^It, up to the pressure where 
helium solidifies Pmax ~ 25atm. Below we present the structure of the vortex based on this EOS, in 
order to give the reader a rough estimate of how realistic the assumption p =const, which we have 
used in this work, is. Because the fluid at T = can only exist in the superfluid phase for pressures 
between SVP and Pmax the density is only allowed to be in the range [po, 1.188po]- From the EOS 
we compute the enthalpy density, w{p) = which for T = can also be regarded as a function 

of p, and is a strictly increasing function. The vortex profile, i.e. p = p{r), is obtained from the 
Bernoulli equation: 

w{p{r)) = w{p^) - 1^^, with v{r) = ^. (93) 

The density (and thus enthalpy) decrease towards the center, and a boundary of the fluid develops at 
r = r^,, with p(r^,) = po- In case we neglected the non-linear terms in the EOS (i.e. had we assumed 



a constant speed of sound), we would have obtained p(r) = pooexp 



, where we recognize 



the dimensionless combination |(— )^ in the exponent. Let us call such a case the simplified model. 



Generally, the speed of sound decreases towards the core of the vortex, and therefore the ratio 



2^ r 

c(r) 

grows monotonically. It is of interest to determine the highest possible value of this parameter; a 
value greater than 1 would signify the presence of an ergoregion. Let poo = xpo; in the simplified 
model we find 

r* = -7^- (94) 
V ^ log X 
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FIG. 11: The profile of density (solid) and speed of sound (dashed) for a vortex in superfluid helium-4 
modeled with the equation of state of [39[. The chosen pressure at infinity is the maximal allowed pressure, 
and the vortex carries the single circulation quantum. Density is shown w.r.t. the saturated-vapor-pressure 
density, psvp = Po- Distances are given in [cm], and the innermost part should not be taken literally, 
because it is of the size of the van der Waals radius for ^He. 

While this result is valid for any fluid with the simplified equation of state p{p) = c^^p (cqo a 
constant), it is curious, that the maximal available x for superfluid helium is 1.188, making r^, ^ l.Tr^. 
Therefore, even though the size of the vortex, r^,, can be increased to avoid the difficulty with the 
van der Waals radius by increasing the circulation k (that is: by trapping more vortices on the 
central wire), the ratio f(r)/c(r) will never exceed 1/1.7 ~ 0.59. No higher value of the ratio can 
be produced with superfluid helium not contained by external walls, within the simplified model 
(including the dependence of c on p does not change this conclusion significantly). 

If, on the other hand, the fluid were contained by an external cylinder of radius R (at pressure 
^ 25atm), it would be possible to have almost arbitrarily large v/c. The values r* ^ Tc are achieved, 
however, only for relatively small diameters R; for instance in order to get r^, = one would need 



to take 



56] X = 1.65 leading to R = 2.359 x 10 ^cm and r^, = 2.175 x 10 ^cm. Even though both of 



these radii can be appropriately scaled up by choosing a larger circulation k, R will always remain 
only by 8% larger than r^,. 

The saturated- vapor-pressure density pq, as well as the maximal density at which the superfluid 
remains liquid, xpo, are clearly thermodynamic properties of the substance under consideration. It 
is difficult to give a principal reason why their ratio should not lead to 1/a/2 logx < 1, although this 
is not the case in superfluid helium-4. 
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Appendix A: Spectral properties of selfadjoint operators in Krein spaces 



Let if be a selfadjoint operator on a Krein space with the indefinite product (., .) (we drop the 



index "H here). The following properties hold (see |27|,|40|,|41|): 

a) if cJ 7^ uj' then the corresponding eigenvectors are orthogonal, (\l/a;, ^oj') = 0, 

b) for each eigenvalue u there also exists an eigenvector to the eigenvalue oJ, 

c) if an eigenvalue u is not real, then the corresponding eigenvector is of zero norm, {^ui-, '^lo) = 0, 

d) if to a certain eigenvalue u there also exists an associated eigenvector, = uj'^'^ + 
(where is the normal eigenvector to u), then the normal eigenvector is of zero norm, 
(^,,^,) = 0. 



According to [27|], in the case considered there (electrostatic potential ^(x)), the eigenfunctions 
together with the associated eigenfunctions form a complete set, so that any element of the phase 
space can be decomposed as (see eq. (3.33) of |28|) 

$(t) = J2 {a^^..e-"^'' + k [K, - e-^'^'*} , with = *^(x) (Al) 



Appendix B: Positivity of the energy functional 

We introduce the decomposition gab = —NaNb + e^ef + ■yab and Ta = aNa + /3e*, (for vortex 
spacetime /3 = and a = 1), and examine the integrand of the energy-functional E[ip]. Using the 
notation: = N'^da, 9$ = e%da, (d^i')'^ = '-)°'^dai>dbi>., we find that 

T^bN'-T' = I [{d^^y + 2£(9^^)(a$^) + {d^tpy + [d.tljf] , (Bl) 

which is positive only as long as \(3/a\ < 1. As this condition is not fulfilled in the ergoregion (in 
our case: for r < 1), the density which is integrated over to obtain the energy can be made 
negative for some classical field configurations. 

Appendix C: Relation to general theory of Bloch/Floquet/Mathieu 

The equation 

^Dl-Ml + K'il + j,)]i; = (CI) 
upon substitution u = \iax, u E [— oo, oo] becomes the standard modified Mathieu equation 

ijjuu - [M^ - cosh(2M)]^ = 0, (C2) 
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which in turn results from the Mathieu equation 



i),, + [M^ - cos(2z)]V' = 



(C3) 



taken along the pure imaginary axis z = iu. In the literature 2J] the following definitions of 
constants in the Mathieu equations are employed: 



= a, K'^ = q. 



(C4) 



This is especially convenient for the reason, that the bound state problem (which appears in case 
there stands a —1 instead of +1 in the regular bracket of the original equation) is just the original 
problem taken along the line(s) u = w ± ij, and with the parameter q = ^iK"^ (the x— equation 
fulfilled by the bound state is recovered by substitution x = e"'). Our study of solutions of the radial 
equation in the space of essential parameters is thus seen equivalent to the study of the solutions of 
the Mathieu equation: 

• for real and positive a's and g's (corresponding to positive and - scattering states for 
non-absorbed rays), for z = iu 

• for real a's and g's with a negative a (corresponding to = — /x^ < - scattering of "absorbed 
rays") z = iu 

• for real a and purely imaginary q (bound states) with z = iu^i^. The wavefunctions need to 
vanish at a; — )■ +cxd, which corresponds to \ni{z) +oo. 

bound state domain 



-ijt/4 



regular Mathieu domain 



ijt/4 



X^modified Mathieu domain 

FIG. 12: Relation between Mathieu functions, and solutions of the radial equation for real frequencies uj. 
The region z — >• ioo, corresponding to x — t- oo, is where the bound state functions should exponentially 
vanish. 

The double eigenvalue problems associated with the regular or modified Mathieu equation possess 
a number of special distinguished solutions. In particular, there is a well-known family of even/odd 
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solutions of period tt or 27r (they are the Mathieu functions). These solutions are defined by Fourier 
series, e.g. 

ce(z) = c„ cos(nz), (C5) 
where one immediately establishes a recursion relation for the c— coefficients, 

V2n C2n = C2n+2 + C2„-2, = [o + {2nf]/q, (C6) 

which lead to a continuous-fraction relation for the ratios = C2n+2/c2n, together with the first 
ratio Gq = a/q. In turn this establishes a family of curves in the (a, g) which are allowed for this 
type of (even-, vr-periodic) solution. Besides the periodic solutions, there are also the so-familiar 
Bloch solutions (in solid state these correspond to "electrons moving through the crystal"), 

^(z) = e'P' cos{nz), (C7) 

where we have used the tilde to emphasize, that a new set of c„ needs to be obtained for each 
quasi-momentum p. The solution to the double eigenvalue problem (family of curves in the (a-q) 
plane) is provided from two sets of continued fractions, Gm = Cm+2/cm, Hm = Cm-2/cm fulfilling the 
condition GqHq = 1. 



A further large family of almost distinguished solutions is provided (see (24|, §20.4) by modi- 
fications of the Floquet (Bloch) solutions, with the same set of coefficients Cm, multiplying the 
functions 

M^) = (-x)^ Besselfc(2v^x), x = [^^^^^l • (C8) 

[cos(2; + 0) 

The coefficient b in the independent variable, as well as the particular type of the Bessel function, 
may be chosen freely. These new types of functions, with different behavior in the complex z plane, 
has, due to the recursion relations, exactly the same spectrum in the (a-q) plane. 
Still further solutions, employing the same coefficients Cm are provided by 

^ 5^52n(-)"BesseU,+,(v^e^^)J„_,(v^e-*^). (C9) 



C2s ^7 

with an arbitrarily chosen s, arbitrary Bessel function Bessel, and the Floquet (Bloch) momentum 
P- 

If, for a given boundary condition, some of the solutions of the above type can be shown to have 
the right behavior at the boundary, than the whole spectrum (family of lines in {K"^, M^) plane) 
is identical to the corresponding Floquet (Bloch) spectrum (with the index p). This would be so, 
because the spectrum is determined by the coefficients c„, and they, in turn, depend only on p (and 
the type of Floquet solution used). However, at this point none of these standard solutions of the 
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Mathieu equation appears to play a role for the vortex-bound-state problem. On the other hand 
standard approximative methods associated with the theory of Mathieu equations 2S, |42| lead to 
efficient approximations for the positions of bound states and for the scattering characteristics, in 
particular in the case of large M^. 



Appendix D: Note about the employed numerical method 

The Mathieu-nature of the sole non-trivial differential equation of this paper makes the use of 
numerical methods unnecessary. However, it proved more efficient to actually use them when dealing 
with the equation (due to the boundary conditions in a form not standard for Mathieu problems, and 
due to the robustness of numerical methods for solving ODEs). The radial equation has therefore 
been implemented in Mathematica, and solutions for arbitrary uj,p,m were obtained starting from 
the appropriate boundary condition at small r. Positions of bound states were obtained by searching 
for zeroes of the Wronskian of the numerical solution with the asymptotic solution which is damped 
for r — )■ oo (functions of Bessel-K type, Km{Kx), see I VI A Til . 
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[45] This step is relatively easy for free fields, which is the case we consider here. 
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where is either the solution of the homogeneous equation, or the solution at the previous position 

in the Jordan chain, as the problem has the form {H — uj)ipa+i = ipa- 
[47] Misprints in eq. (26) of [s^ do not influence this result; the power of r multiplying the cosine function 

is special to Case's form of the problem. 
[48] Restoring the arbitrary value of E, the acoustic rays with E'^ ~ would reach the core for J G [—2E, 0] . 
[49] And, of course, by {uo^m) — t- (— i^, —m) symmetry, with negative-frequency, positive-norm solutions. 
[50] Friedman's work is related to the astrophysical context, and therefore assumes asymptotic flatness. 

Although the vortex (acoustic) spacetime due to its cylindrical symmetry does not conform to this 

assumption, we just use it to make the point, that bound states would provide a way out of the 

conclusion of the theorem of Friedman. 
[51] We also extend the domain of m to real numbers. 

[52] An investigation, comparable to the acoustic-spacetime description, would need to use the extended 
f„.al.n, of [22], WhO. th.s, to the be. of out fatowledge, has not been attempted yet. it wonld 
provide an interesting extension of the work of Flaig and Fischer, |13l ]. 

[53] oo = 5 • lO^^cm is the Bohr radius. 

[54] We denote the mass of the boson constituent of the condensate by m^. 
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momenta, as the giant vortices begins to emerge, there remains a fraction of the Tkachenko lattice for 
r > with a constant lattice spacing d. (see 9| ) ; the velocity field corresponding to the missing lattice 
(for r < Rfi) still circulates around the core of the giant vortex. 

[56] Without the exterior container, this value of x would require poo = 212atm (using unjustified extrap- 
olation of the equation of state) . 
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